General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA CONTRACTOR REPORT 


NASACR-141425 

(NA3A-CR- 141 425) INIIIAL CEVhLCPaENI Qi A N77-30443 

MEriKC OF SIGNIFICANT WAVEHtlGHT ESTIKATICN 
FOR b EOS -11 1 (Applied Science Assoc id tes, 

Inc-, Apex, N.C,) 13 p HC A02/HF ACl UncidS 

CSCL 08C 03/3> 42UG3 

INITIAL DEVELOPMENT OF A METHOD OF SIGNIFICANT WAVEHEIGHT 
ESTIMATION FORGEOS-III 


G. S. Hayne 


Prepared Under Contract No. NAS6-2520 by 


Applied Science Associates, Inc. 
105 East Chatham Street 
Apex, North Carolina 27502 


N/\SA 

National Aeronautics and 
Space Administration 

Wallops Flight Center 

Wallops Island Vi'^ginia 23337 
AC 804 824 3411 


r? r 


RECt: 

sp,sA sn FACiua 
^{^PUT BRANCH , 


August 1977 



1. Repoft No. 


2. Government Accession No. 


3. Recipient's Catalog No. 


NASA CR-141425 


4 . Title and Subtitle 

Initial Development of a Method of Significant Waveheight 
Estimation for GEOS-III 


5. Report Date 
August 1977 


6. Performing Organization Code 


7. Author(s) 

G.S. Hayne 

9. Performing Organization Name and Address 

Applied Science Associates, Inc. 
105 East Chatham Street 
Apex, North Carolina 27502 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Wallops Flight Center 
Wallops Island, VA 23337 


8. Performing Organization Report No. 


10. Work Unit No 


11. Contract or Grant No 

NAS6-2520 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

Details are presented for the numerical method by which ocean surface significant waveheight 
estimates are produced from the GEOS-3 radar altimeter data. Four parameters characterizing 
the expected radar mean return waveform are determined from the altimeter's sixteen sample-and- 
hold waveform samplers through use of an iterative, least-squares approach. One of the four 
parameters is a risetime term which has contributions from both the transmitted radar pulse 
width and the rms ocean surface height; the estimated significant waveheight is extracted from 
this risetime. The noise character of the estimation process leads to occasional negative 
significant waveheight estimates, and the origin of these non-physical negative results is 
discussed. Possible modifications and areas for additional investigation are indicated. 


17. Key Words (Suggested by Author(s)) 

18 Distribution Statement 

GEOS-3 

Unclassified - unlimited 

sea state 
radar altimeter 

STAR Category 35 

■ 1 — ■ ■ 


19. Security Oeuif. (of thi* report) 20 Security Classif (of this page) 21 No of Pages 22. Pnte 

Unclassified Unclassified 11 


For sale by the Netionel Techni-al (nformation Service, Springfield. Virginia 22151 















TABLE OF CONTENTS 


Page 

1.0 INTRODUCTION AND SUMMARY OF RESULTS 1 

2.0 SIGNIFICANT WAVEHEIGHT ESTIMATION 1 

2.1 Background 1 

2.2 Development of Estimation Method 2 

3.0 CONCLUSIONS, RECOMMENDATIONS, AND FUTURE PLANS 9 

REFERENCES 10 


ii 



1.0 INTRODUCTION AND SUMMARY OF RESULTS 


Work accomplished under this portion of Contract NAS6-2520 comprises 
derivation of an approximate numerical technique for use in estimating sig- 
nificant waveheight from the GEOS-III radar altimeter waveform data. 

The expression developed for been incorporated into the GEOS- 

III data processing system at Wallops Fllgh.: Center, and was utilized in a 
near-realtime mode during the February 1976 North Atlantic High Sea State 
Exercise just concluded. 

The algorithm currently in use represents only the result of an Initial 
development designed for early computerization so that the data processing 
system would be complete. Work on refinement of this method as well as 
other approaches to estimation is still in progress under follow-on 

tasks. 

2.0 SIGNIFICANT WAVEHEIGHT ESTIMATION 

2.1 Background 

The feasibility of estimating sea surface significant waveheight 
using short pulse radar altimetry has previously been established [1]. Theo 
retical limits on termination resolution have also been established 

[ 2 ]. 

The algorithm developed in this task for estimation from GEOS-III 

waveform data is an extension of the analytical work previously done [3]. 
This prior work established a means of describing an altimeter average re- 
turn waveform using the four parameters: A, t^, 3t and C* where 

A was a general waveform amplitude factor, 

t was a reference time, 
o 

3 was a pulsewldth-rlsetlme parameter, and 

^ was antenna off-nadir pointing angle. 

The parameter 6 Included effects of both the altimeter transmitted 
pulse width and the rms ocean surface roughness. The effect of changes in 
ocean surface roughness Is primarily to alter the slope of the average 
ocean-scattered radar return's leading edge; the present estimation 
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algorithm has been developed from this emphasis on slope of the rise-time 
portion of the return. For the GEOS-III altimeter, the pointing angle ef- 
fects in the leading edge of the return can be ignored because of the rela- 
tively broad beamwidth of the antenna pattern; in this case the more com- 
plicated mean return waveform of Reference 3 and other work, reduces to simply 
an integrated Gaussian (the Probability Integral). One additional parameter 
has been added, a factor to account for a possible non-zero baseline in the 
altimeter's waveform sampler-recorder average return signal. The prob- 

lem has then been reduced to one of determining four parameters characteriz- 
ing the mean return waveform using suitable time-averaging of the 16 wave- 
form S&H gates on GEOS-III, followed by recovery of the sea-state-caused 
portion of the risetime parameter. 

Since changes in surface roughness alter the slope of the average re- 
turn leading edge, only the rise-time portion of the return need be consid- 
ered in estimating surface roughness. The algorithm developed for this es- 
timation, therefore, is derived on this basis. Thus it involves the four 
parameters a, b, c, and d, where: 

a = return amplitude 
b = time origin 
c = return rise-time, and 
d * return baseline. 

For the GEOS-III altimeter, the pointing angle effects in the leading edge 
of the return can generally be ignored due to the relatively broad beamwidth 
of the radar antenna pattern. It is necessary to introduce a factor which 
accounts for a possible nonzero baseline in the sampler recorded average re- 
turn - this is the parameter "d". 

2 . 2 Development of ^ ^ Estimation Method 

Based on an iterative, least squares approach [4], a means of estimat- 
ing the parameter "c" (from which can be derived) is presented in the 

following; equations 1-6 and accompanying discussion are for the general 
4-parameter least-squares problem. 

Then, let there be a four-parameter theoretical function, y^, which is 
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a function of the Independent variable t; 


- y(a,b,c,d,t) 


( 1 ) 


Write the first-order Taylor series approximation to this (at a point 
a^,b^,c^, and d^ in the parameter space) as follows; 


V i'b® * ''J >’d* 


( 2 ) 


where 


y » y(a ,b ,c ,t) 

•'O "^000 


^a da 


a j b f 
o o 

c ,d , 
o o 


db 


a y b f 
o o 

o o 
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'c dc 


o o 

C »d„.t 
o o 


dd 


o o 

c ,d ,t 
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and a ■ a-a , 3 * b-b , y ■ c-c and 6 ■ d-d , Note that the factors y . , 
o o ' o o 

l«a, ...,d, are functions of t. 
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Then form a sum over all N values of the errors squared; 




( 3 ) 


Minimizing E with respect to a , 


^ - 0 - 2 '£ 


W 


or 


«Ey3 “ ^^a^yr^o^ 


( 5 ) 


Similar equations follow from ^ ^ ^ , and all four such equations 

de dy d6 


rewritten in matrix form give; 
2 
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^Vc 


^''ayb 

r 2 

^>'b>'c 

^>'b’'d 


=^>'b>'c 

^>-0 



“\ 
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6, 


Ey^Cyj-yo) 
2yb(yi-y„) 
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( 6 ) 


Returning now to the specific estimation algorithm for which the 

specific is taken as the Probability Integral on a baseline, let 


't ' ^ “* 


where P(') is the Probability Integral, P(z) -y Z(q)dq, 

-oo 


where Z(q) is the Gaussian function, Z(q) 


- q^/2 


/2tt 


and the parameters a,b,c,d, have the physical meanings 



a “ amplitude 
b ■ time origin 
c • rlsetlme 
d “ baseline 

In terns of Z, P as defined above. 



and y. “ 1 

a 


Rewriting the 4-matrlx from the least-squares non-linear problem, 



Substituting in the values of y tyvty ly. etc., 

abed 
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'(B[- v(^) - v(^)- -.VV, . 

i*:l •- -* i“l 
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Completion of the problem requires finding the inverse of the 4 by 4 

matrix on the left side of (8) above in order to solve for a, B, y. and 5. 

The initial guesses a , b , c , and d are corrected to produced revised 
® o o o o 

estimates a, b, c, d since a*a +a, b*b +B, c = c +>, and d = d +6. 

o o o o 

The sign is used to emphasize that a a; since only the first order 

terms were retained in the original Taylor series at equation (2), and a 
by definition was a = a^+6, thus a is only an approximate estimate of a. 
Similarly b is an estimate of b, c of c, and 3 of d. 


In this iterative procedure, the next step then is to use a, b, c, d 
in place of the origlual a^, b^, c^, d^, and repeat the entire cycle to 
produce (presumably) improved estimates of a, b, c, d. The elements through 
•^4 above are of course re-evaluated each time through the cycle; the indi- 
cated sums are for the N*16 separate values from the 16 GEOS-III waveform 
samplers. The iterative procedure is repeated unti* a defined relative 
error criterion is satisfied (junless the matrix inversion fails first and 
causes an error exit); a typical loop exit criterion is to require that the 
total sum of errors squared decrease by less than 0.1% of the sum of errors 
squared in the last previous time through the loop. For a reasonable first 

guess a ,b ,c ,d (and "reasonable" first guesses can be as much as a factor 
o o o o 

of two away from the final answer in our experience with the particular 
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functional form in this ^-related problem), the problem will converge in 
typically two to four iterations to reasonable final answers for a,b,c, and d. 
Given the risetime parameter, c, determined by this iterative least-squares 
procedure above, the estimation of follows: 

The convolutional model for the mean sea-scattered pulse return assumes 
that the risetime parameter c is the composite result from the calm sea (as- 
sumed Gaussian) equivalent radar pulsewidth ("equivalent" implies inclu- 
sion of trcinsmltter and receiver bandwldths, etc.) and the rough-sea rms 
heights o Because of the convolution of Gausslans of widths 0 and a , 

2$ SC 

c^ » (9) 

s c 


or 


a 

s 




Here O is two-way ranging time in units of nanoseconds. Conversion of 
s 

ranging times (nanoseconds) to surface height (meters) is accomplished by 
multiplying by the two-way speed of light “ .15 m/ns. Since is four 

times the rms surface height relative to mean sea level, must also be 
multiplied by four; 

= 4(.15) Vc^ - (10) 

V ~2 2 

c - in m, c and in ns). 

Expression (10) is true only in a limiting sense. Actually the c in 
(10) is an experimental estimate based on a finite data span, and as such, 
c is noisy and fluctuating. The quantity under the radical in (10) is con- 
sequently capable of going negative in some data spans and thus ending the 
problem since any practical computer program will terminate upon trying to 
take the square root of a negative number. To avoid this, (10) is modified to 

Hi/ 3 - 0.6^ |c^ -0^1 ^ Sign(Ic^-oJ^ I) . (11) 


*Tracking loop jitter has been ignored. 
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Our bringing out the sign from within the radical was done in order 

2 

to allow recovery of the quantity c (knowing already the value of O ) 

c 

for any possible c value o to The actual value of used in expres- 
sion (11) was based upon early observations of several passes over known 
calm seas; the value currently used is * 8.55 ns, but this should be 
reexamined now that considerably more data are available. 

3.0 CONCLUSIONS, RECOMMENDATIONS, AMD FUTURE PLANS 

The algorithm just described has been utilized in the GEOS-III 

data processing programs at Wallops Flight Center; specifically it has been 
installed in the computer program GAP more fully described in Leitao ^ al . 
[5]. Work now underway is comparing this product with other ground 

truth efforts particularly for the February 1976 North Atlantic data. 

Additional work is now being done on the algorithm itself in two 

areas. First, the results can be smoothed over a several data frame 

period by either averaging the individual-frame results or by averag- 

ing together the waveform sample -and -ho Id gate results over the several- 
frame interval and using these to produce a single value representing 

the several-frame period. Averaging the single-frame is 

simpler to do but we are now examining carefully the question of whether 
this estimate converges to that from the wavefoi i averaging approach. 

Second, the effects of tracker jitter are now being examined; these 
were left out of the abvive-described first version of a procedure. 

The outcome of this additional work will be reported later. 
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